On locations and properties of the multicritical point of Gaussian and ±J Ising spin 
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We use transfer-matrix and finite-size scaling methods to investigate the location and properties 
of the multicritical point of two-dimensional Ising spin glasses on square, triangular and honeycomb 
lattices, with both binary and Gaussian disorder distributions. For square and triangular lattices 
with binary disorder, the estimated position of the multicritical point is in numerical agreement with 
recent conjectures regarding its exact location. For the remaining four cases, our results indicate 
disagreement with the respective versions of the conjecture, though by very small amounts, never ex- 
ceeding 0.2%. Our results for: (i) the correlation-length exponent i' governing the ferro-paramagnetic 
transition; (ii) the critical domain-wall energy amplitude rj; (iii) the conformal anomaly c; (iv) the 
finite-size susceptibility exponent 'j/i'; and (v) the set of multifractal exponents {rik} associated to 
the moments of the probability distribution of spin-spin correlation functions at the multicritical 
point, are consistent with universality as regards lattice structure and disorder distribution, and in 
good agreement with existing estimates. 

PACS numbers: 75.50.Lk, 05.50.+q 
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I. INTRODUCTION 

In this paper we investigate quenched random-bond 
Ising spin - 1/2 models on regular two-dimensional lat- 
tices, namely square [SQ], triangular [T], and honeycomb 
[HC]. For suitably low concentrations of antiferromag- 
netic bonds, it is known that such systems exhibit ferro- 
magnetic order at low temperatures. We consider only 
nearest-neighbor couplings Jij, with strengths extracted 
from identical, independent probability distribution func- 
tions (PDFs). We specialize to the following two forms 
for the latter: 



P(Jy)=P<5(J,, 
P(Jy) = J- exp 



27rcr 



Jo) + (1 ~ p) S{J^j + Jq) 

{Jjj - Jq) 
2cr2 



(Gaussian) . (1) 



Our units are such that Jq = 1 in the former case, and 
(7 = 1 in the latter. A critical line on the T — p (±J), 
or T — Jo (Gaussian), plane separates paramagnetic and 
ferromagnetic phases; a spin-glass phase for comparable 
amounts of plus and minus couplings is absent here, on 
account of the systems under consideration being two- 
dimensional. 

For general space dimensionality d> 2 there is a sec- 
ond line of interest on the temperature-disorder plane, 
along which the internal energy has a simple analytic 
expression, and several exact results have been derived, 
known as the Nishimori line (NL)^. The shape of the 
NL is known exactly, and given by 



-2/T 



1-p 



T = — 
Jo 



{p>l: 



(Gaussian) . 



(2) 



The intersection of the ferro-paramagnetic boundary 
with the NL is a multicritical point^., the Nishimori point 



(NP). A conjecture regarding the possibly exact location 
of the NP has been put forward, which invokes the effects 
of duality and gauge symmetry arguments on the repli- 
cated partition function of quenched random Zq mod- 
els4i.fi.i2i. With further extensions to consider non self-dual 
latticesii^, numerically exact predictions have been pro- 
duced, for the Z2 (Ising) model, for all lattices and inter- 
action distributions considered here. Versions of the con- 
jecture adapted for hierarchical lattices have been con- 
sidered as well^. 

Locations of the NP predicted by the conjecture gener- 
ally agree very well with results obtained by other means. 
However, the remaining discrepancies provide compelling 
evidence that, at least in some cases, the conjecture may 
not be exact. First, on the SQ lattice, several very ac- 
curate numerical estimates for the ± J coupling distribu- 
tion place the conjectured location^, Pc — 0.889972 . . . , 
outside the corresponding error bars (though it differs 
from the central value typic ally by less than 0.1%). One 
has: Pc = 0.8906(2) (Refs. 0.8907(2) (Ref. HI), 

0.89081(7) (Ref. II4). For a Gaussian distribution, the 
conjecture givea^ifi. Jqc = 1.021770, while Ref. [ill finds 
Joc = 1.02098(4). Second, it has been shown that the 
exact renormalization-group solution for three pairs of 
mutually dual hierarchical lattices disagrees with the per- 
tinent form of the conjecture, by up to 2%^. 

Very recently, these issues have been addressed via the 
proposal of an improved conjecture, first applied to hier- 
archical latticesii, and later extended to regular ones^^. 
Broadly, this corresponds to considering duality prop- 
erties appHed to a (usually small) cluster of sites on the 
lattice under examinatio n^^i^^ , as opposed to the original 
conjecture which deals only with the partition function 
of a single bond (the principal Boltzmann factor )4. The 
improved conjecture predicts the location of the NP to 
be well within the error bars of recent numerical work 
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for the SQ ±J case: an average over four slightly dif- 
fering implementations gives Pc = 0.89079(6), though so 
far disagreement persists for the Gaussian distribution, 
as the improved estimate is Jqc — 1.021564iS.. For hi- 
erarchical lattices, the gap between conjecture and ex- 
act renormaHzation-group solutions has essentially been 
bridged by the new approach- 

Existing numerical results for T and HC lattices (±J 
distribution only)^'^^ broadly agree with an early form 
of the original conjecture, applicable to pairs of dual lat- 
tices^: with the binary entropy 

H{p) = -p\og^p~ {I ~ p) \og^{l - p) , (3) 

it is predicted that, for a pair of mutually-dual lattices 1 
and 2, 

Hi2 = Hipic) + H{p2c) = 1 . (4) 

Ref.0 finds Pc = 0.835(5) and 0.930(5), respectively for 
T and HC, which implies 0.981 < H12 < 1-042; these 
estimates were refined in Ref. m to Pc = 0.8355(5) [T] 
and 0.9325(5) [HC], giving H12 = 1.002(3). 

Further developments^ii^ enabled the production of 
pairs of individual predictions (always obeying Eq. 
with a suitably-adapted form of Eq. ^ for the Gaus- 
sian case). In the framework of the original conjecture, 
these are: Pc = 0.835806 [T] and 0.932704 [HC] (±J)^; 
Joc = 0.798174 [T] and 1.270615 [HC] (Gaussian)!^. For 
the improved conjecture (±J only), two slightly differing 
implementations give the pairs: pc — 0.835956 [T] and 
0.932611 [HC]; pc = 0.835985 [T] and 0.932593 [HCp8. 

Here, we numerically estimate the location and critical 
properties of the NP on the T and HC lattices. For the 
± J case, we refine the results given in Ref. [l^, checking 
our data against the more stringent predictions of Refs.0, 
[15I : for the Gaussian distribution, we are not aware of any 
existing results, apart from those given in Ref. [Tsl for the 
conjectured location of the NP. For completeness, and to 
provide consistency checks of our methods, we revisit the 
SQ lattice problem, investigating both distributions. 

We apply numerical transfer-matrix (TM) methods 
to the spin-1/2 random-bond Ising model, on strips of 
SQ, T, and HC lattices, of widths 4 < iV < 14 sites 
(SQ), 4 < iV < 13 sites (T) and 4 < TV < 14 sites 
(even values only, HC). We take long strips, usually 
of length M = 2 X 10^ columns (pairs of columns for 
HC, because two iterations of the TM are needed to 
restore periodicity). For each of the quantities evalu- 
ated here, averages (Q) are taken over, and fluctuations 
(AQ)i.ms calculated among, Ng independently-generated 
samples, each of length M. As discussed extensively else- 
whereH, the sample-to-sample fluctuations (AQ)rms vary 
with M~^/^, and are essentially TVs-independent, pro- 
vided that Ns is not very small. The averaged values 
(Q) themselves still fluctuate slightly upon varying Ns, 
but the corresponding fluctuations A(Q) die down with 
increasing Ng. We found that, for M as above, mak- 
ing Ns = 10 already gives A(Q)/(AQ)„„s < 0.1, thus 



Table I: Intervals Ap, AJo scanned along the NL in our cal- 
culations, for lattices and coupling distributions [binary (± J) 
and Gaussian (G)] specified in column 1. Afp gives total num- 
ber of pairs {p,N) or {Jo,N) at which quantities of interest 
were calculated. See text. 



Type Ap, AJo TVp 

SQ, ±J 0.8868 - 0.8948 123 

SQ, G 1.00125 - 1.03875 140 

T, ±J 0.830 - 0.842 126 

HC, ±J 0.9266 - 0.9386 86 

T, G 0.7794- 0.8169 134 

HC, G 1.254 - 1.287 94 



this constitutes an adequate compromise between accu- 
racy and CPU time expense. Typical upper bounds for 
(AQ)rms/(Q) are lO"** for free energies, 10"'^ for domain- 
wall energies (see Section lUl below for deflnitions) . 

We scanned suitable intervals of p or Jq along the NL, 
spanning conjectured and (when available) numerically- 
calculated positions of the NP, as shown in Table HI For 
a given lattice and interaction distribution, we took sam- 
ples at Np — Np{N) equally-spaced positions for each 
lattice width A^, generally starting with Np > 18 for 
small A^, and decreasing to iVp = 9 for TV > 8, giving 
the totals denoted by Afp = J2n ^pi^) in Table H 

The Mersenne Twister random-number generator- 
was used in all calculations described below. In all cal- 
culations pertaining to the ±J disorder distribution, a 
canonical ensemble was used, i.e. for a given nominal 
concentration p of positive bonds, these were drawn from 
a reservoir initially containing aipNM units (a^ = 2,3,3 
respectively for i = SQ,T, HC). This way, one ensures 
that fluctuations in calculated quantities are consider- 
ably smaller than if a grand-canonical implementation 
were use d^^'^^ . 

In Sec. ini domain-wall energies are computed, and 
their flnite-size scaling allows us to estimate both the lo- 
cation of the NP along the NL, and the correlation- length 
index, yt = 1/v which governs the spread of ferromag- 
netic correlations upon crossing the ferro-paramagnetic 
phase boundary. The conformal anomaly, or central 
charge, is evaluated in Sec. IIIII In Sec. ITVl uniform sus- 
ceptibilites are calculated, and the associated exponent 
ratio ^y/iy is evaluated (for Gaussian coupling distribu- 
tions only). In Sec. [Vl we specialize to T and HC lat- 
tices, with Gaussian disorder distributions, and investi- 
gate the moments of assorted orders of the probability 
distributions of spin-spin correlation functions. Finally, 
in Sec. IVIj concluding remarks are made. 



II. DOMAIN- WALL SCALING 

For a strip of width L, in lattice parameter units, of a 
two-dimensional spin system, the domain- wall free energy 
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(Tl is the free energy per unit length, in units of T, of a 
seam along the full length of the strip. For Ising spins, 
(7l = fL~ fb^ '^ith f£ (ff) being the corresponding free 
energy for a strip with periodic (antiperiodic) boundary 
conditions across. Within a TM description of disordered 
systems, ctl = — ln(AQ /A^) where InA^, InAg are the 
largest Lyapunov exponents of the TM, respectively with 
periodic and antiperiodic boundary conditions across. 

The duality between correlation length ^ and interface 
tension a is well-estabHshed22 for pure two-dimensional 
systems, and carries over to disordered cases. In a finite- 
size scaling (FSS) contexli^i, this means that ctl must 
scale with 1/L at criticality, a fact which has been used in 
previous studies of disordered systems^^, including inves- 
tigations of the NP— iiiii^. Prom conformal invariance^^ 
one has, at the critical point: 

LaL ^TTT] , (5) 

where, for pure systems, 77 is the same exponent which 
characterizes the decay of spin-spin correlations. In the 
presence of disorder, however, the scaling indices of the 
disorder correlator (i.e., the interfacial tension) differ 
from those of its dual, the order correlator (namely, spin- 
spin correlations}^. Nevertheless, the constraints of con- 
formal invariance still hold, thus the amplitude of the 
domain wall energy remains a bona fide universal quan- 
tity^i. Por the NP, recent estimates on the SQ lattice 
(±J couplings) give ?? = 0.691 (2)i2J^i2i.. 

We have calculated A^, A^ for strips of SQ, T 
and HC lattices. Recalling that both L in Eq. ^ 
and the correlation length ^ (of which the surface 
tension is the dual) are actual physical distances 
in lattice parameter units^i^>22ii^, one finds (see 
Ref. [iB) that, in terms of the number of sites N 
across the strip, the appropriate expressions for the 
scaled domain-wall energy are of the form: ryjv = 
m{T,z) = Q{N/TT){{lnAP{T,z)-\nA^{T,z)), with 
Q = 1, 2/V3, \/3/2 respectively for i = SQ, T, HC, and 
where z — p (± J) or Jq (Gaussian). At (Tc, Zc) one must 
have limAT^oo Vn — ^7, the latter being a universal quan- 
tity. 

Close to the multicritical NP, the scaling directions 
are respectively the NL itself, and the temperature 
^-jj.jg3.i3,29 Therefore (neglecting corrections to scaling), 
along the NL the single relevant variable corresponds to 

Z - Zc. 

According to finite-size scaling^!, the curves of scaled 
domain-wall energy calculated for different values of 
N, T, z along the NL should then coincide when plotted 
against x = N^^'^ (2 — Zc). 

Bearing in mind that corrections to scaling may be 
presen t ^^1^^ , we allow for their effect from the start. 
Thus, we writei^: 

m = f[N'/'' (z - Zc)] + N-'^g[N'/'' {z ~ z,)] , (6) 

where w > is the exponent associated to the leading 
irrelevant operator. Close enough to the NP the scal- 
ing functions in Eq. ([6]) should be amenable to Taylor 



expansions. One has: 

j = l k=0 

(7) 

We adjusted our TM data to Eq. by means of mul- 
tiparametric nonlinear least-squares fits. The goodness 
of fit was measured by the (weighted) per degree of 
freedom (x^^ ^ f ). We tested several assumptions on km, 
jm, oj, via their effect on: (i) the resulting X^d o f ' (ii) 
the stability of the final estimates for z^, rj, l/i^, and (iii) 
the broad compatibility of estimates for 77 and l/v with 
existing results for assorted two-dimensional lattices and 
coupling distributions (under the reasonable assumption 
of universality, which is however provisional, and must 
be weighted against the bulk of available evidence) . 
We found that: 

(1) a parabolic form, jm = 2, is adequate for the descrip- 
tion of the broad features of data, similarly to conclusions 
from the Monte-Carlo study of Ref. [13; 

(2) Neglecting corrections to scaling (all bk = 0) generally 
gave a X^d.o.f . least one order of magnitude larger than 
if such corrections are incorporated; 

(3) Fixing km = and allowing uj to vary gave a final es- 
timate uj ^ 0.1 — 0.2, which is too low to qualify as a bona 
fide correction-to-scaling exponent; the same happens if 
one allows fc^ > 1 with a variable uj; 

(4) For fixed w, using km = ^ reduces the x^d o f by 
a factor of 2-3 compared with making km = 0, while 
no noticeable improvement is forthcoming from allowing 
/cm > 1, again in line with Ref. [T3: 

(5) For fixed uj between 1 and 2, results for r) and 1/v 
are in fair accord with point (iii) above; ; also, for this 
range of w, x^d.o.f. is minimized, at ~ 0.1 — 0.2, compared 
to any alternative combination of fixed and variable pa- 
rameters described in this paragraph. The coexistence of 
these facts indicates that, within the assumed scenario 
of describing corrections to scaling via a single (effective) 
exponent, the range of uj just quoted is the one that op- 
timizes a universality-consistent picture. 

Thus, we kept jm = 2, km = 1, allowing 1 < < 2 in 
what follows. Results for uj = 1.5 are shown in Table Ull 

Since the error bars quoted in the Table only refiect 
uncertainties intrinsic to the fitting procedure, we now 
illustrate (see Table Hill below) the quantitative effects of 
relaxing some of the assumptions specified above. This 
is especially important as regards Zc, whose calculated 
fractional uncertainty is one to two orders of magnitude 
smaller than those for 1/z/, 77. Additional checks on the 
robustness of such narrow error bars are therefore in or- 
der. 

For instance, considering the T, ± J case, fixing w = 
1, 2 gives respectively pc = 0.83611(8), 0.83565(6), with 
X^d o f varying by less than 10% against its value for 
UJ — 1.5. Overall, it seems that a realistic error bar 
should at least include the fitted values of Zc obtained for 
UJ = 1 and 2. Table ITTTl shows such estimates, denoted by 
zj^™, where the associated uncertainties refiect the spread 
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Table II: TM estimates of critical quantities Zc (2 = p, Jo), 
1/v, and rj for lattices and coupling distributions [binary (±J) 
and Gaussian (G)] specified in column 1. Column 2 gives 
conjectured values of «c; quotations from Refs.QQQ, and (for 
T and HC lattices with G distribution) Ref.lTsI refer to original 
conjecture, while all others refer to improved conjecture. All 
fits used u! = 1.5 (fixed), see Eq. ^ and text. 



Type 


conj. 


Pc, Joc 




V 


2 

X d.o.f. 


SQ, ±J 


0.8899731"^ 


0.89061(6) 


0.64(2) 


0.689(2) 


15/116 




0.89079(6)" 










SQ, G 


1.021770^"^ 


1.0193(3) 


0.65(3) 


0.680(2) 


28/133 




1.021564i^ 










T, ±J 


0.835806^ 


0.83583(6) 


0.65(2) 


0.691(2) 


18/119 




0.83597(2)'' 










HC, ±J 


0.932704** 


0.93297(5) 


0.65(1) 


0.702(2) 


15.5/79 




0.93260(1)'' 










T, G 


0.798174i^ 


0.7971(2) 


0.66(2) 


0.689(2) 


17/127 


HC, G 


1.270615i^ 


1.2689(3) 


0.64(3) 


0.690(2) 


11/87 



"average over four values from improved conjecture, Ref. llE 
''average over two values from improved conjecture. Refill 



Table HI: For lattices and coupling distributions [binary (± J) 
and Gaussian (G)] specified in column 1, column 2 gives crit- 
ical quantities Zc"'' {z = p, Jo), averaged over values from fits 
for uj = 1 and 2 (see text). Coefficients 60, 61 (see Eq. ^) fit- 
ted for u) = 1.5; the index (0) for the last two columns denotes 
quantities obtained from fits where corrections to scaling were 
neglected. 



Type zr bo 

SQ, ±J 0.89065(20) -0.126(3) 

SQ, G 1.0193(4) 0.009(3) 

T, ±J 0.83588(23) -0.145(3) 

HC, ±J 0.93300(15) -0.142(3) 

T, G 0.7972(6) -0.106(3) 

HC, G 1.2691(10) -0.152(4) 



bi 



.(0) 



2(0) 
A d.o.f. 



-3.7(4) 0.8898(1) 416/118 
-0.35(18) 1.0195(1) 30/135 
-1.4(4) 0.8348(1) 520/121 
-2.5(4) 0.9322(1) 499/81 
-0.31(18) 0.7948(3) 163/129 
-0.64(15) 1.2635(9) 325/89 



between these extreme values (their own intrinsic uncer- 
tainties generally being somewhat smaller, see above and 
Table [TT| . A remarkable exception is the SQ, G case, 
for which the estimate of Joc is virtually unchanged as 
UJ varies in the range described. This instance is also an 
exception in that the amplitude of the correction term 
60 (column 3 of the Table) is much smaller than for all 
other cases; consequently, neither Jqc nor the (resp. 
columns 4 and 5) change appreciably when corrections to 
scaling are ignored. The latter is not true for any of the 
other cases studied. 

Our assessment of the estimates quoted in Table [llll for 
the location of the NP is as follows. 

For SQ, ±J our results are in agreement with the 
improved c onject ure^, and with numerical data from 



improved coniecturei*, 
Refs. Hliillilil. For 



T, ± J our range of estimates is 



roughly consistent with the conjecture, both in its orig- 
inal^ and improvedi^ versions. It is also at the upper 
limit of the early estimate Pc = 0.8355(5)^^. 

For all remaining cases, our numerical data indicate 
that the conjecture fails to hold, albeit by rather small 
amounts, 0.2% at most. For the ±J distribution, our 
results for both SQ and HC indicate that the conjectured 
position of the NP lies in the paramagnetic phase (for SQ, 
this is true only for the original conjecture). On the other 
hand, for the Gaussian distribution and all three lattices, 
according to our estimates the conjecture places the NP 
slightly inside the ordered phase. 

For HC, ± J the result in Table [TO is again at the 
upper end of the range given in Ref. UM, Pc = 0.9325(5). 
Note also that our estimate for SQ, G lies farther from 
the conjecture than the numerical value given in Ref. [ill, 
namely Jqc = 1.02098(4) (thus, this latter also places 
the conjectured location of the NP inside the ordered 
phase). The above estimates of Pc and Jqc for T and HC 
lattices, when plugged into Eq. (|4]), using Eq. ^ and its 
counterpart for Gaussian distributional^, result in: 

i7(pi,) + H{p2c) = 0.9986(12) (± J) ; 
F(Joci) + H{Jqc2) = 1.0014(10) (Gaussian) , (8) 

both narrowly missing the conjecture of Eq. (jj]). 

As regards the correlation-length exponent and the 
critical amplitude 77, we found that, for each lattice and 
coupling distribution, the error bars quoted in Table 
are wide enough to accommodate the variations in cen- 
tral estimates, both when one sweeps uj between 1 and 
2 as above, and when Zc is varied between the Hmits es- 
tablished in Table IHII No evidence emerges from the 
data which justifies challenging our earlier assumption 
of universality. From unweighted averages over the re- 
spective columns of Table Ull we quote = 1.53(4), 
r] = 0.690(6). These are to be compared to the recent 
results f = 1.50(3) (SQ, Ref. [13), 1.48(3) (SQ, Ref.[ii|), 
1.49(2) (T and HC, Ref.[lB), 1.527(35) (SQ, Ref.[l2Lall 
for ±J distributions; see also 1^ = 1.50(3) (SQ, Ref [Till. 
Gaussian. For the critical amplitude, we recall (all for 
±J): Tj = 0.691(2)i2^^ (SQ); 0.674(11) (T), 0.678(15) 
(HC), both from Ref . [H . 

The overall quality of our scaling plots is illustrated 
in Figures [T] and [2l We chose to display data for T and 
HC, Gaussian distributions, because for these there are 
fewer data available in the literature. As the last column 
of Table [ll] shows, the Xd oi remains very much in the 
same neighborhood for all cases studied. 



III. CENTRAL CHARGE 

We used the free-energy data generated in Section |ll] 
also to estimate the conformal anomaly, or central charge 
c, at the NP. This is evaluated via the finite-size scaling 
of the free energy on a strip with periodic boundary con- 
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Figure 1: (Color online) Triangular lattice, Gaussian coupling 
distribution: scaling plot of domain-wall free energies, sub- 
tracting corrections to scaling: 77^" = rjN — N~" g(N^^'^ {Jo — 
Joc)) [see Eq. ((GJ ], against the scaling variable N^^'^ (Jo — Joc). 
Error bars are smaller than symbol sizes. The full line is a 
quadratic fit to corrected data. 
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Figure 2: (Color online) Honeycomb lattice, Gaussian cou- 
pling distribution: scaling plot of domain-wall free ener- 
gies, subtracting corrections to scaling: v'n'^'^ = Vn — 
N~" g{N^^'^ {Jo — Joc)) [see Eq. ((S} ], against the scaling vari- 
able N^^'^{Jo — Joc). Error bars are smaller than symbol sizes. 
The full line is a quadratic fit to corrected data. 



Table IV: Conformal anomaly c and non-universal higher- 
order coefficent d, from fits of critical free-energy data to 
Eq. ((QJ. Last column gives fitted values of c under the ass- 
sumption that d = in Eq. 



Type 


c 


d 


c [d = 0] 


SQ, ±J 


0.463(3) 


0.13(1) 


0.478(2) 


SQ, G 


0.461(4) 


0.14(3) 


0.476(2) 


T, ±J 


0.459(3) 


0.01(1) 


0.461(1) 


HC, ±J 


0.457(5) 


0.02(2) 


0.462(2) 


T, G 


0.454(4) 


0.06(3) 


0.461(1) 


HC, G 


0.468(15) 


-0.05(6) 


0.459(5) 



ditions acrosi 



,30 



/(T„7V) = /(T„oo) + ^ + A + o 



1 



(9) 



where /(Tc, 00) = liniL^oo f{Tc,L) is a regular term 
which corresponds to the bulk system free energy. For 
disordered systems, Eq. ^ is expected to hold when 
the configurationally-averaged free energy is considered, 
with c taking the meaning of an effective conformal 
anomalj*ii>^>2^ . 

By writing only even powers of iV~^ in Eq. |[9]), it is 
assumed that only analytic corrections come up^. While 
this is true, e.g., for pure Ising systems, a counterexam- 
ple is the three-state Potts ferromagnet for which free- 
energy corrections in TV"^"", N~^'^° . . . , ujq = 4/5, are 
presen t^^i^^ . Although not much is known about the op- 
erator structure at the NP, existing central charge es- 
timates in this case have been derived via Eq. ^ so 
far with fairly consistent results, namely c=0.464(4) ^°i^^ , 
0.46(1)25.. We shall return to this point at the end of this 
Section. 

We have evaluated free energies at the predicted loca- 
tions of the NP given in Table IIIH both at the central 
estimates and at either end of the respective error bars. 
We found that such values can be calculated with suf- 
ficient accuracy, via interpolation of those already com- 
puted at the sets of equally-spaced points used originally 
in Section ini Results for the central charge are displayed 
in Table ITVl where error bars for all quantities mostly re- 
flect uncertainties intrinsic to the fitting procedure itself, 
as our estimates are rather stable along the predicted in- 
tervals of location of the NP. Indeed, it is expectedii that 
at criticality the calculated conformal anomaly passes 
through a maximum as a function of position along the 
NL. 

In line with earlier findings^^, one sees that for SQ and 
both coupling distributions, ignoring the fourth-order 
term in Eq. ^ shifts the final estimate of c by some 
4 — 5 error bars, away from the expected universal value 
~ 0.46. On the other hand, for T and HC the fitted d is 
much closer to zero than for SQ; furthermore, for these 
latter lattices, results obtained fixing d = appear gen- 
erallly more consistent with universality, and with less 
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Table V: For the zero-field susceptibility and lattices as spec- 
ified in column 1 (all with Gaussian coupling distributions), 
columns 2, 3 , 4 give the leading correction-to scaling expo- 
nent as fitted, and the corresponding ^jv and X^d.o.f. ; columns 
5, 6 give the two latter quantities, now taken by keeping ijj 
fixed during the fitting procedure, and averaging over result- 
ing values for a; = 1 and 2 (see text). 



Type 






X^d.o.f. 




^ ,2 avc 
A d.o.f. 


SQ, G 


1.3(3) 


1.79(2) 


0.127 


1.793(6) 


0.129 


T, G 


0.7(3) 


1.81(1) 


0.118 


1.814(6) 


0.128 


HC, G 


0.4(6) 


1.79(4) 


0.17 


1.804(7) 


0.18 



spread, than those found with d kept as a free parame- 
ter. Overall, we interpret the above results as indicating 
that: (i) there is no evidence for universality breakdown 
as regards the conformal anomaly; taking this as true, 

(ii) there appears to be no unusual (non-analytic) free- 
energy scaling correction N~'^° with 2 < ijJq < 4; and 

(iii) it is possible that the fourth-order term is d = for 
T and HC, similarly to the case of pure Ising systems^^. 

IV. UNIFORM SUSCEPTIBILITIES 

We calculated uniform zero-field susceptibilities along 
the NL, for SQ, T and HC lattices and only for Gaus- 
sian distributions, as done in previous investigations for 
_[_ ji6,36 ^ For the finite differences used in numerical dif- 
ferentiation, we used a field step 6h = 10~^ in units of 
T. We swept the same respective intervals of Jq quoted 
in Table m 

Finite-size scaling arguments^! suggest a form 

XN^N'^/'^fiN'/'^iJo-Joc)], (10) 

where xn is the finite-size susceptibility, and 7 is the 
susceptibility exponent. In order to reduce the number 
of fitting parameters, we kept l/v and Jqc fixed at their 
central estimates obtained in Sec. [ill and allowed j/f to 
vary. Again, we took corrections to scaling into account. 
Following Ref. [l^, we write: 

In X = J In TV + ^ a^. {J„ - J^^y N'/- + 

km 

+N-^Y.^k{Jo-hc)^N''''' . (11) 

Similarly to Section [Ul above, we found that choosing 
jm = 2, fc„i = 1 enables one to obtain good fits to nu- 
merical data, with xl o f — ^l- ~ 0.2. The consequences 
of keeping w as a free parameter or, on the other hand, 
fixing its value during the fitting procedure, can be seen 
in Table [Vl While the fitted value of lu for SQ looks ac- 
ceptable, the same cannot be said of that for HC, as the 



associated error bar allows even slightly negative values 
(the result for T being half-way between the other two). 
Also, by keeping w as a free parameter, one gets an error 
bar for ^/v that is at least twice that obtained if to is 
kept fixed between 1 and 2, without any noticeable im- 
provement in the x^d o f ■ the other hand, using fixed 
oj above this latter range results in a slow but steady loss 
of quality: for example, for the T lattice, lj = A gives 
X^d of— Thus, although the idea of allowing uj to 
vary freely seems, in principle, the correct thing to do, 
the results in this particular case do not appear to be 
obviously more reliable than those averaged for fixed uj 
between 1 and 2. We then decided to use these latter 
as our main reference. Taking an unweighted average 
over the three estimates for (7/1^)'*™ gives the final value 
^/v = 1.804(16). This is to be compared to the following 
(all for ±J distributions): 1.80(2)|^ [SQ]; 1.795(20) [T] 
and 1.80(4) [HC], both from Ref. [H; 1.80-1.82" [SQ]. 

V. CORRELATION FUNCTIONS 

Our study of correlation functions is based on previous 
work for SQ21, T and HC latticesi^ (±J only). In this 
Section, we specialize to the Gaussian distribution, for T 
and HC lattices only. On the NL, the moments of the 
PDF for the correlation function between Ising spins en, 
(Tj are equal two by tw o^i^i^^i^^ : 

[(a,a,)2^-i] = [(a,a,-)2^] , (12) 

where angled brackets indicate thermal average, square 
brackets stand for configurational averages over disorder, 
and £ = 1, 2, ... . 

At the NP, conformal invarianceiS is expected to 
hold, provided suitable averages over disorder are consid- 
eredi2^>i^>2i^did2d3 . On a strip of width L of a square 
lattice, with periodic boundary conditions across, the 
disorder-averaged fc-th moment of the correlation func- 
tion PDF between spins located respectively at the origin 
and at (x, y) behaves at criticality as: 

^ z""" , z = [sm\?{TTx/L)+fiin^{'Ky/L)Y'^ . 

(13) 

For the T and HC lattices, the same is true, provided 
that the actual, i.e., geometric site coordinates along the 
strip are used in Eq. (fTSl . Details are given in Ref. [l^. 
Note that Eq. I|12p implies 772^-1 — ^21- 

As in earlier worki^^, we concentrate on short- 
distance correlations, i.e., where the argument z is 
strongly infiuenced by y. Such a setup is especially conve- 
nient in order to probe the angular dependence predicted 
in Eq. lfT3|) . which constitutes a rather stringent test of 
conformal invarian ce proper ties. 

Following Refs. Ild|la|37l . we extract the decay-of- 
correlations exponents rjk, via least-squares fits of our 
data to the form nik ^ 2;"'"=. We also consider the ex- 
ponent ryo which characterizes the zeroth-order moment 
of the correlation-function distribution's., i.e. it gives the 
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typical, or most probable, value of this quantity (see, e.g., 
Ref. I43I and references therein). One has, in the bulk. 



Table VI: Estimates of exponents rj k , from least-squares fits of 
averaged moments of correlation-function distributions to the 



Go{R) = cxp [ln(croO'i?,) 



R 



-Vo 



(14) 



Calculations on strips of the ± J SQ lattice, at the early 
conjectured location of the NP— , gave the estimate rjo = 
0.194(1)^ . 

As seen in earlier work-^, for strip widths = 10 or 
thereabouts, finite- width effects are already mostly sub- 
sumed in the explicit L (i.e., N) dependence of Eq. (fT3|) . 
However, some detectable (albeit tiny) variations in the 
calculated values of averaged moments of the correlation 
function PDF may still persist upon varying N. These 
are of course minimized at the critical point where the 
bulk correlation length diverges. We have calculated cor- 
relation functions for < 12, for values of Jo within the 
error bars given for the location of the NP in Table IIIII . 
We have seen that along these intervals of Jo, the trend 
followed by such averaged moments against A^- variation 
is as follows: for T, it cannot be distinguished from stabil- 
ity within error bars, while for HC it is slightly downward 
(of the order of one error bar from A^ = 10 to A' = 12). 
For fixed N and Jq, one error bar associated to intrinsic 
fiuctuations is < 1%. 

Table IVTl gives numerical results of the fits for fc = 0, 
and odd A: > 1 (we have also calculated even moments 
for k > 2 and checked that Eq. (fT2|) holds). One sees 
that T and HC estimates are quite consistent with each 
other for all k. On the other hand, for fc = 0, 1, and 3 
they fall slightly below their existing counterparts, given 
in columns 4, 5, and 6 of the Table. For k — 5 and 7, as 
a consequence of generally wider error bars, all estimates 
are broadly compatible with one another. 

Physically, obtaining (via least-squares fits) a smaller 
[larger] than expected value for the decay-of-correlations 
exponent would indicate that it is being evaluated inside 
the ordered [paramagnetic] phase, instead of right at the 
critical point. 

Applying these ideas to the fc = case, we recall that 
the result of Ref. [H for SQ, ± J was calculated at the 
originally conjectured position of the NP— . By now, it 
seems well established that this point is in the disordered 
phase (see Table [ill. Therefore, the value of Ref. [H 
should be taken as an upper bound, which is obeyed by 
our present estimates. 

For fc = 1 and 3, one might use the same argument as 
above to argue that the result of Ref. [sl is too large, as 
it was calculated at the same point as that of Ref. Is^. 
On the other hand, this cannot be said of the addi- 
tional estimates quoted in the Table, all of which are also 
larger than ours (though in some cases the respective er- 
ror bars overlap, or at least touch each other). Using 
the reasoning described above, one would infer that for 
T and HC with Gaussian distribution, the ranges of lo- 
cations for the NP given in Table IHII are in fact both 
inside the ordered phase. Since these latter, in their 
turn, put the conjectured NP position also inside the 
ordered phase, the final conclusion would be that the ac- 



form nik 



for z defined in Eq. I|13p . For columns 2 and 



3 (this work), central estimates and error bars reflect averages 
between results for A = 10 and 12, as well as variations from 
scanning Jo along the error bars for locations of NP given in 
Table Hill Columns 4, 5, 6 quote existing data for comparison. 
For SQ, all results are for ±J coupling distribution, unless 
otherwise noted. 

k T (G) HC (G) T {±J)^ HC {±J)^ SQ 

0.185(3) 0.184(3) - - 0.194(1)3^ 

1 0.178(2) 0.178(2) 0.181(1) 0.181(1) 0.1854(17)^ 

0.1854(19)1-2 

0.183(3)1^ 
0.1848(3)ii 
0.1818(2) [G)ii 
0.180(5)i^ 

3 0.250(2) 0.252(2) 0.251(1) 0.252(1) 0.2556(20)^^ 

0.2561(26>iS 
0.253(3)i^ 
0.2552(9)ii 
0.2559(2) [G}ii 

5 0.296(2) 0.300(5) 0.297(2) 0.296(2) 0.300(2)^ 

0.3015(30)ifl 
0.3004(13)ii 
0.3041(2) [G)ii 

7 0.331(4) 0.336(6) 0.330(2) 0.329(3) 0.334(3)^ 

0.3354(34)ifl 
0.3341(16)ii 
0.3402(2) [G)ii 



tual location of the NP differs from the conjecture by an 
amount larger than predicted by domain- wall (DW) scal- 
ing: Jo'^^^i < Jo°w < The slight downward trend 
against increasing A^, reported above for HC, would be 
consistent with this scenario. However, we have not seen 
a similar trend for T. 

One should note also that all the discrepancies re- 
marked upon are rather small: the single worst 
regards central estimates, is that of the present result 
ryi = 0.178 against 771 = 0.1854-ft t ^T^ amounting to 4%, 
or ~ 3.5 times the respective error bar. Given that the 
quoted values (especially those for the associated uncer- 
tainties) are likely to depend on details of the respective 
fitting procedures, the resulting picture looks mixed. 

In conclusion, existing evidence does not seem strong 
enough to state that our estimates from Sec. Ull for the 
location of the NP on T and HC lattices (Gaussian dis- 
tribution) are definitely inside the ordered phase. 
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VI. DISCUSSION AND CONCLUSIONS 

We have used domain- wall scaling techniques in SecUll 
to determine the location of the Nishimori point of Ising 
spin glasses. In the analysis of our data we allowed for 
the existence of corrections to scaling, see Eqs. ([6]) and 
0. 

Results for the SQ lattice, ± J distribution, show that 
such corrections play a crucial role in the finite-size scal- 
ing of domain-wall energies. Indeed, when they are 
taken into account, the estimated position of the NP is 
Pc — 0.89065(20), in excelent agreement with recent and 
very accurate numerical worl i^^i^^ . One can see from the 
two last columns of Table lTTll that. if corrections to scaHng 
are ignored, the value of Pc which minimizes the Xd o f 
(though at a level ^ 30 times that obtained when cor- 
rections are incorporated) is instead Pc = 0.8898(1), very 
close to the original conjecture and incompatible with 
the above-mentioned body of numerical evidence. In ret- 
rospect, one sees that the domain- wall scaling result of 
Ref. [H for this case, namely pc — 0.8900(5), essentially 
suffers from the effect of ignoring corrections to scaHng 
(though even so it still picks out the correct exponent, 
l/y= 1.45(8)ift). 

Going over to SQ, Gaussian, domain-wall scaHng for 
strips of widths 4 < < 14 sites gives Joc = 1.0193(4), 
lower than both the conjecture (original and improved) 
and the result of Ref. [O, -he = 1.02098(4)). In that 
Reference, the mapping of the spin problem to a network 
model (described in Ref. [13) enabled the authors to reach 
significantly larger lattice sizes than here. The result just 
quoted was obtained by extrapolation of 11 < iV < 24 
crossing-point data, without explicit account of correc- 
tions to scaling (which, as those authors show, do pro- 
duce a trend reversal around N = 8, and are expected 
to have negligible effect for the large widths used in the 
extrapolation) . It may be that our own data fail to incor- 
porate an underlying trend which only comes about for 
larger systems. Nevertheless, the stability of our results 
for this particular case is remarkable, as pointed out in 
the initial discussion of Table IIIII 

Our results for T and HC, ±J distribution, are 
marginally compatible with, but more accurate than, the 
earlier ones of Ref. \1Q, though for T they are also broadly 
compatible with the conjecture in both original and im- 
proved versions, for HC our estimate in Table ITTll lies at 
least two error bars away from the conjecture. 

For T and HC, Gaussian distribution, in both cases 
the discrepancy between our results and the conjecture 
is again of the order of two error bars. 

Consequently, as shown in Eq. we predict the 
duality-based conjecture of Eq. |[4]) to be narrowly missed, 
for both ± J and Gaussian cases, though on opposite sides 
of the hypothesized equality. 

As regards the exponent ly and critical amplitude t] (see 
Eq. ^) which are also estimated via domain- waH scaling, 
we have found no evidence of nonuniversal, lattice- or 
disorder distribution- dependent, behavior. Therefore, 



from unweighted averages over all six cases studied, we 
quote I' = 1.53(4), rj = 0.690(6). Both are in very good 
agreement with existing numerical results (see the end of 
Sec. mi for detailed comparisons). 

The conformal anomaly values calculated in Sec. ITTTl 
are in good agreement among themselves and with pre- 
vious estimatesiSiiii'^. Our fits for the non-universal co- 
efficient d of the fourth-order correction to the critical 
free energy suggest that d = for T and HC lattices 
(while definitely d for SQ). This would be similar 
to the lattice-dependent structure of corrections for pure 
Ising systema^^. An unweighted average of values from 
Table ITVl (using results of fits with d ^ for SQ, and 
with d = for T and HC) gives c = 0.461(5). 

In Sec.[lV]we evaluated uniform zero-field susceptibil- 
ities, by direct numerical differentiation of the free en- 
ergy against external field. Only Gaussian distributions 
were considered, for SQ, T, and HC. Though our results 
show some lattice-dependent spread, the error bars for 
{l still overlap in pairs. It is known that suscepti- 
bility calculations are prone to larger fiuctuations than, 
e.g., domain- waH energy oneaii. In the Monte Carlo sim- 
ulations of Ref. this effect was reduced by considering 
the quantity (where ^ is the finite-size correlation 

length), which behaves more smoothly than x on its own. 
Our final estimate (averaged over results for all three 
lattices), ^/v = 1.804(16), compares favorably (albeit 
somewhat close to marginally) with the corresponding 
one from Ref.0, -f/v = 1.820(5). 

Finally, in Sec. |V] we applied conformal-invariance 
concepts to the statistics of spin-spin correlation func- 
tions, extracting the associated multifractal scaHng ex- 
ponentsi^'iiii^'^i. We only examined T and HC lattices, 
for Gaussian coupling distributions. The overall picture 
summarized in Table [Vll points towards universality of 
the exponents {r]k}, though some smaH discrepancies re- 
main. The case fc = 1 is especially relevant, on ac- 
count of its connection with the uniform susceptibility 
via the scaling relation 7/1/ = 2 — 771. While our re- 
sult ryi = 0.178(2) is somewhat lower than existing data 
from direct calculations of correlation functions, it gives 
^jv = 1.822(2) when inserted in the scaling relation. 
This agrees very well with the above-quoted estimate^^, 
7/z/= 1.820(5). 

In summary, we have produced estimates of the loca- 
tion of the NP on SQ, T and HC lattices, and for ±J 
and Gaussian coupling distributions. Though these are 
consistent with existing conjectures for SQ and T (both 
±J), they appear to exclude the respective conjectured 
values for the remaining cases. However, the discrepan- 
cies are very small, amounting to 0.2% in the worst case 
(SQ, Gaussian). Furthermore, we have assessed several 
critical quantities (amplitudes and exponents) , and found 
an overall picture consistent with universality as regards 
lattice structure and disorder distribution. 
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